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ABSTRACT 


The  dynamics  of  the  Free  Electron  Laser  are  governed  by  Maxwell’s  equations  which  causes 
many  highly  nonlinear  regimes  to  exist  in  Free  Electron  Laser  Physics.  This  thesis  will  examine  two 
such  areas  and  develop  simple  models  to  describe  the  highly  dynamic  and  rich  behavior  two  of  these 
regimes. 

In  the  strong-field,  high-current  regime,  the  Free  Electron  Laser  driving  current  can  be 
modeled  by  a  single  macroparticle  representing  the  trapped  electrons.  When  the  trapped  electrons  act 
collectively  as  a  macroparticle,  solutions  which  include  synchrotron  oscillations  can  be  found  for  the 
self-consistent  pendulum  and  wave  equations. 

In  an  FEL  oscillator  with  low  single-pass  gain,  the  evolution  of  the  optical  wave  can  lead  to 
sideband  development.  This  phenomenon  is  studied  by  applying  Maxwell’s  equations  to  an  oscillator 
with  two  optical  modes  and  deriving  a  two-mode  wave  and  pendulum  equation.  The  two-mode  wave 
and  pendulum  equations  are  implemented  numerically  on  computers  so  that  the  onset  of  the  sideband 
can  be  explored. 
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I  INTRODUCTION 


Free  Electron  Lasers  (FELs)  are  proving  to  be  an  exciting  and  practical  source  of  radiation 
with  a  wide  assortment  of  applications  over  a  broad  range  of  the  RF  spectrum.  First  proposed  in 
1970  by  John  Madey  [Ref.l],  and  experimentally  demonstrated  by  Madey  and  his  colleagues  in  1975 
and  1976  [Ref.  2,3],  they  are  now  being  considered  for  use  in  a  wide  variety  of  applications.  The 
FEL  is  based  on  the  original  work  by  Motz  who  demonstrated  that  light  could  be  amplified  by  an 
electron  beam  propagating  through  a  periodic  magnetic  field  [Ref.  4,5,6].  In  1960,  the  ubitron  was 
developed  by  Phillips  which  operates  on  similar  principals  to  the  FEL  [Ref.  7]. 

The  FEL  has  many  attributes  that  make  it  attractive  for  a  wide  assortment  of  uses.  The  most 
notable  of  these  attributes  is  the  FEL’s  continuous  tunability.  FELs  have  been  demonstrated  to 
operate  over  a  range  of  frequencies  that  vari^  by  a  factor  of  10  [Ref.  8].  Closely  related  to 
continuous  tunability  is  the  "designability"  of  the  FEL.  The  same  basic  theory  and  mechanisms  that 
work  for  a  FEL  operating  in  the  millimeter  range  will  also  work  in  the  tens  of  nanometer  range. 
FELs  can  also  be  scaled  to  high  power  applications  since  the  gain  medium  contains  only  a  relativistic 
electron  beam  in  a  periodic  magnetic  field,  and  does  not  suffer  the  nonlinear  material  affects  that 
plague  convention  atomic  lasers  in  the  high  power  regimes.  This  also  makes  FELs  very  reliable  since 
the  weakest  link  in  FEL  design,  the  accelerator,  is  generally  considered  very  reliable.  FELs  do  suffer 
two  major  drawbacks,  size  and  cost.  Some  of  the  ciurent  work  in  FELs  is  aimed  at  reducing  the  size 
which  will  help  with  the  associated  high  cost  [Ref.  8]. 

Another  active  area  of  FEL  research  is  in  the  high-power  regime  where  they  have  many 
applications  ranging  from  plasma  heating  to  weaporuy.  Unfortunately,  in  the  high-power  regimes. 
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two  different  areas  where  nortlinearides  arise,  and  develop  simple  theories  that  provide  useful  insight 
into  these  dynamic  regimes. 

Chapter  II  will  provide  the  necessary  background  on  classical  FEL  theory.  Chapter  IQ  will 
develop  a  simple  theory  that  treats  the  bunched  electrons  in  the  FEL  as  a  single  macroparticle.  This 
greatly  simplifies  the  FEL  dynamics  and  allows  simple  solutions  to  be  derived  describing  FEL 
behavior  in  this  highly  nonlinear  regime.  The  solutions  can  then  be  used  to  study  the  trapped-particle 
instability,  sideband  development,  and  chaos  in  the  FEL.  The  solutions  can  also  be  used  for  design 
criteria.  Chapter  IV  develops  a  simple  two-mode  model  of  the  FEL  oscillator  by  solving  Maxwell’s 
equations  for  two  optical  modes  in  the  oscillator.  The  model  is  then  used  to  study  the  growth  of 
sidebands  in  an  oscillator.  A  quantitative  understanding  of  this  phenomenon  will  add  a  tremendous 
degree  of  control  in  FEL  oscillator  design. 
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II  FREE  ELECTRON  LASER  THEORY 


A.  BASIC  FREE  ELECTRON  LASER  PHYSICS 

An  FEL  consists  of  two  major  components,  an  electron  accelerator  to  produce  a  relativistic 
electron  beam  and  a  wiggler  or  undulator  with  a  static,  periodic  magnetic  field  to  "wiggle"  the 
electrons  with  a  periodic  oscillatory  motion.  It  is  this  periodic  motion  of  the  electrons  in  the 
undulator  that  causes  radiation  in  the  forward  direction.  Some  of  this  spontaneous  radiation  is  saved 
in  the  laser  resonator  causing  stimulated  emission  that  leads  to  coherent  radiation.  Since  the  electron 
beam  is  highly  relativistic,  the  electrons  "see"  the  undulator’s  magnetic  field  wavelength  as  being 
much  shorter  due  to  relativistic  Lorentz  contraction.  Subsequently,  the  optical  field  produced  by  the 
undulating  electrons  is  of  a  much  shorter  wavelength  than  the  undulator’s  magnetic  field  wavelength. 

The  above  is  a  Quantum  Mechanical  description  of  the  basic  FEL  physics.  Despite  this,  the 
FEL  mechanism  is  essentially  classical  and  can  be  adequately  described  with  Electro-Magnetic 
Theory.  Figure  2-1,  illustrates  the  major  featiures  of  the  classical  approach. 

The  incoming  electrons  enter  the  undulator  and  begin  to  oscillate  due  to  the  periodic  magnetic 
field  of  the  undulator  as  shown  in  the  top  of  Figure  2-1.  The  oscillating  action  of  the  electrons  causes 
them  to  emit  photons  in  a  manner  analogous  to  a  radiating  antenna.  In  the  middle  of  Figure  2-1,  one 
undulator  wave  section  is  enlarged  to  show  both  tlie  optical  and  magnetic  fields  of  the  undulator 
superimposed  over  the  oscillating  path  of  the  electron.  The  forces  acting  on  the  electrons  are 
governed  by  the  Lorentz  force  law 

>  0  X  ,  y  =  — ,  2-1 

at  me  me 

where  ^  =  vie  and  y  =  1  -  0^ .  The  speed  of  light  is  c ,  the  mass  of  the  electron  is  m ,  the  electron 
charge  is  e ,  and  the  velocity  of  the  electron  is  i^.  £  and  S  represent  the  electric  and  magnetic  fields 
present  in  the  undulator  from  both  the  optical  wave  and  the  static  undulator  magnetic  field.  By 
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Figure  2-1:  Major  components  of  the  Free  Electron  Laser  illustrating  the  fundamental  physics  of  the 
laser. 
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examining  the  bottom  third  picture  in  Figure  2-1,  it  can  be  seen  that  the  electrons  are  effected  by  the 
optical  Held  and  the  effect  is  dependent  on  the  electron’s  phase  relative  to  the  optical  field.  This 
coupling  between  the  electrons  and  the  optical  wave  is  best  illustrated  with  an  example.  When  the 
electron  has  a  velocity  component  both  in  the  z  direction  and  out  of  the  page,  and  the  magnetic  field 
of  the  optical  wave  is  directed  downward,  the  electron  will  experience  a  force  retarding  its  forward 
motion  and  causing  the  electron  to  lose  energy.  Conservation  of  energy  requires  that  the  electron’s 
lost  energy  go  to  the  optical  wave  causing  subsequent  amplification  of  the  optical  field.  Different 
electron  phases  with  respect  to  the  optical  field  can  also  cause  the  electron  to  gain  energy  from  the 
optical  field.  Therefore,  careful  analysis  is  required  of  the  electron  and  optical  field  coupling  to 
produce  amplification  of  the  optical  field. 

B.  ELFXTRON  DYNAMICS  AND  THE  PENDULUM  EQUATION 

In  order  to  fuUy  understand  the  interaction  of  the  electrons  and  the  optical  field  inside  the 
undulator,  the  equations  governing  the  evolution  of  the  electrons  must  be  developed.  Begin  by 
assuming  a  helically  polarized  static  magnetic  undulator  field  of  the  form: 

B  =  B(cos(^ot),  sin(^),  0)  ,  2-2 

where  B  is  the  field  strength,  z  is  the  distance  alon^  he  longitudinal  axis  of  the  undulator,  ond^g 
is  the  wave  number  of  the  undulator’s  magnetic  field.  Also  assume  that  a  circularly  polarized  plane 
wave  is  present  in  the  undulator  cavity  of  the  form 

Eg  -  £(cosi|f, -sintj;,  0),  Bg  =  £(sini)i,  cosi)i,0)  ,  2-3 

where  E  is  the  optical  field  strength  and  =  fe  -  tor  +  (|)  with  <o  =  Jtc .  The  wave  number  for 
the  optical  field  is  denoted  by  k,  time  by  t ,  and  the  optical  phase  by  <|> . 
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Substituting  equations  2-2  and  2-3  into  the  Lorentz  force  equations,  2-1,  and  separating  the 


electron  motion  into  their  transverse  and  longitudinal  components  the  following  is  obtained 

ti(Yp  )  « 

— — ^  =  -~[£(l-pp(cos\}f.-shn|i,0)+p^(-sinA^,cos^,0)l  , 
at  me 


2-4 


and 


dt 


-_L[£(P  cosy  -  P  siny)  +  5(P  sin/rpZ  -  P  cosItoZ)]  , 


2-5 


_  =  -_f-£ip  cos\|/-P  siny]  ,  2-6 

dt  me  '  ^ 


where  =  (pj,,Py,0).  Since  the  electrons  are  highly  relativistic,  the  approximation  1-P,  *  0 
can  be  used  to  show  that  £(1  -  p^)  <  flp^  [Ref.  8].  This  simplifies  the  transverse  equation  2-4 
and  allows  it  to  be  solved  by  inspection  giving 


=  -—{co^,sink^,0)  ,  2-7 

where  K  =  eBX^Ilnmc^  is  a  dimensionless  parameter  that  characterizes  the  properties  of  the 
undulator.  In  a  FEL  K  is  usually  close  to  1.  The  transverse  motion  solution,  2-7,  assumes  perfect 
beam  injection  so  that  there  are  no  constants  of  integration.  The  transverse  motion  is  on  a  much 
smaller  scale  than  the  average  z  motion  along  the  z  axis  [Ref.  8]. 

Substituting  the  transverse  motion,  2-7,  into  equation  2-6  gives 


Y 


eKE 

- - (cos(^„r)cosy  -  sin(^„z)siny))  . 

mcf 


2-8 


Using  trigonometric  identity  for  the  cosine  of  the  sum  of  two  angles  and  the  definition 
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2-9 


^  B  {k  *  k^)z  -  (Ot  the  change  of  the  electron’s  energy  is  found  to  be 

Y  =  -f^cos(C  +  <|>)  . 

ymc 

^  is  a  dimensionless  "microscopic"  variable  that  describes  the  electron’s  phase  with  respect  to  the 
optical  wave  and  the  undulators  magnetic  field.  ^  on  the  "scale"  of  X  ,  the  wavelength  of  die 
optical  wave  in  the  undulator,  since  k>k^  and  z  >=  ct.  The  exchange  of  energy  of  the  electrons  is 
easily  understood  in  terms  of  the  electron’s  phase  from  2-9.  For  example,  electrons  with  phases  such 
that,  -n/2  ^  +  <t>)  ^  ti/2  ,  will  receive  energy  from  the  optical  wave  causing  them  to  accelerate. 

Other  electrons  at  phases,  +  «|>):S37t/2  will  lose  energy  to  the  optical  wave.  Since  some 

electrons  gain  and  some  electrons  lose  energy,  a  beam  of  initially  random,  but  essendally 
continuously  spaced  electrons,  will  begin  to  bunch  during  the  trip  down  the  undulator.  This  bunching 
process  is  the  important  feature  of  the  FEL  gain  mechanism  as  will  be  discussed  in  more  detail  in 
section  D. 

To  further  simplify  the  dynamics  of  the  electrons  in  the  undulator,  it  is  convenient  to  relate 
Y  to  ^  .  Using  the  definition 

Y  =  .  1  -  P’  =  1  -  p;  -  Pi  ,  2-10 

and  recalling  2-7,  for  the  motion,  we  have 

Y  ’  =  1  -  p;  -  ,  2-11 

Y 

since  .  Solving  equation  2-11  for  P^  and  differentiating  with  respect  to  time  gives 

p  =  2-12 

using  Y  >  1 ,  and  P^  «  .  Differentiating  the  definition  of  ^  twice  respect  to  time,  and  solving 
for  P^  gives 
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Equating  2-11  and  2-12  gives 


d  =_J_ 

'  (* » 


2-13 


1  .  ^  f"  .  2-14 

Y  (A:  +  kjc  I  + 

When  one  wavelength  of  light  passes  over  an  electron  as  the  electron  traverses  one  undulator 
wavelength,  the  electron  is  in  resonance  with  the  optical  wave,  and  the  electron  is  oscillating  at  the 
same  rate  as  the  optical  field  so  that  t  -  0  ■  Using  this,  the  resonance  condition  can  be  derived 
which  will  be  used  to  simplify  2-14.  Starting  with  the  definition  of  ^  and  differentiating  with  respect 
to  time  yields 

C  =  (A:  +  ^)cp^-fcc  =  0  .  2-15 

Solving  2-11  for  and  substituting  into  2-15  yields 

y'i  2-16 

^  ~  "oTF) 

since  y  >  1  • 

Combining  equations  2-14  and  2-16  and  solving  for  y  gives 


Cy  K  ^  Cy  * 

{k  +  k^c  2X  {k  +  k^c  2k^ 


2-17 


Since  k»k  , 

O  ' 


Y 


-  Ji 

2k^c 


yC 

2u)„ 

0 


2-18 


Equating  equation  2-9,  the  energy  equation,  and  equation  2-18,  and  solving  for  ^  gives 
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2eu>XE 

C  =  — ; - cos(C  +  4)) 

Y  me 


2-19 


which  is  the  simple  pendulum  equation.  To  put  the  pendulum  equation,  2-19,  in  dimensionless  form, 
define  a  dimensionless  time  x  ■  ctjL  ,  such  that  0  ^  x  ^  1,  where  L  is  the  length  of  the 
undulator.  Applying  the  chain  rule  gives  =  t{Llc)^  so  that 

2-20 


2e(oJCEL^  AneKELN 

-cos(g  +  <|))  =  - - — — cos(g  +  <|)) 


y^mc^ 


y^mc^ 


since  kjc  =  oo^  =  c2n/X^  .  iV  is  the  number  of  periods  in  the  magnetic  field  of  the  undulator  and 
the  refers  to  differentiation  with  respect  to  x  ,  as  will  always  be  the  case.  Defining 
\a  1  s  {4neKELN)l(_y^mc^)  as  the  dimensionless  optical  field  strength  so  that  the  optical  field  is 
a  =  |a|c'^  gives 

^  =  lalcos(^  +  «|>)  . 

The  pendulum  equation  given  in  equation  2-21  is  valid  for  both  strong  and  weak-fields.  The 
only  restrictions  placed  on  using  equation  2-21  are  that  y  be  high  enough  that  the  electron  is  not 
bent  excessively  by  the  undulator’ s  magnetic  field  and  that  the  FEL  is  not  so  efficient  that  y 
decreases  substantially  as  the  electron  traverse  the  undulator.  The  optical  field  strength  \a  |  in  front 
of  the  cos  term  explains  why  stronger  optical  fields  produce  rapid  electron  evolution. 

As  the  FEL  normally  saturates,  the  electrons  lose  enough  energy  that  they  drop  out  of  the  gain 
bandwidth.  As  the  electrons  continue  to  lose  energy  and  shift  across  the  gain  bandwidth,  the  beam's 
phase  velocity  v  h  ^  decreases  by  Av  «  AnNAyly  =  -2n  .  By  either  decreasing  the 
undulator ’s  magnetic  field  wavelength  along  the  z  axis,  or  by  decreasing  the  magnetic  field 
strength  B  along  z,  or  both,  the  electrons  could  be  accelerated  in  phase-space  to  maintain  resonance 
with  the  optical  wave.  Each  method  is  conceptually  equivalent  and  both  are  called  "tapering"  the 
undulator  [Ref.  8].  The  phase  acceleration  can  be  added  to  the  pendulum  equation  as  an  added 
constant  that  is  turned  on  near  saturation  at  time  x^.  This  makes  the  pendulum  equation 
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lalcos(C+4>) 

5  +  |a|cos(C  ■*■♦) 


x<x 


2-22 


* 

t>t, 

where  5  is  the  necessary  added  phase  acceleration.  When  the  phase  acceleration  is  created  by 
decreasing  the  undulators  wavelength,  6  -  -2nN^/K^  and  when  the  acceleration  is  provided  by 
weakening  the  magnetic  field  5  ~  -AnNK^6BIB{i*K^)  [Ref.  8].  If  8  >  |a|  ,  then  the  path  of 
the  electron  in  phase-space  will  not  repeat  itself  nor  return  to  its  initial  conditions.  Paths  in  phase- 
pace  that  do  not  repeat  themselves  are  called  open  orbits.  Paths  that  do  close  on  themselves  and 
repeat  are  called  closed  orbits.  If  {n  |  >  8  some  of  the  electrons  can  be  trapped  in  closed  orbits, 
and  with  |a|  >  8  most  of  the  electrons  can  be  trapped  near  resonance  around  ^  n  producing 
high  gain.  Successful  tapering  requires  strong  optical  fields  [Ref.  8]. 

C.  OPTICAL  WAVE  EVOLUTION  AND  THE  WAVE  EQUATION 

In  both  PEL  oscillators  and  PEL  amplifiers,  spontaneous  emission  rapidly  forms  a  classical 
wave  with  some  degree  of  coherence  since  the  bandwidth  of  the  spontaneously  emitted  light  is 
comparable  to  the  inverse  of  the  number  of  undulator  periods  traversed  by  the  electrons  [Ref.  8].  As 
discussed  during  the  pendulum  equation’s  derivation,  there  is  either  energy  loss  or  gain  depending 
on  the  electron’s  phase  relative  to  the  optical  wave. 

The  optical  wave  is  taken  to  vary  slowly  in  time  during  an  optical  period  (£<o>£,  <^<tD(j>) 
and  vary  slowly  in  space  over  an  optical  wavelength  {E'<.kE,  .  This  approximation  is 

known  as  the  slowly  varying  amplitude  and  phase  approximation  (SVAP)  and  is  consistent  with  the 
narrow  bandwidth  expected  for  any  laser.  The  evolving  optical  wave  can  be  viewed  as  a  carrier  wave 
made  up  of  a  single  high  frequency  that  is  modulated  by  a  complex  wave  envelope  that  varies  slowly 
in  amplitude  and  phase  over  many  optical  wavelengths.  It  is  this  complex  envelope  that  meets  the 
SVAP  approximation.  If  the  complex  wave  envelope  did  not  meet  the  SVAP  approximation,  the  light 
from  the  laser  would  have  a  much  broader  bandwidth  than  expected  and  poor  coherence.  This  is 
evident  through  Fourier  analysis. 

To  begin  the  derivation  of  the  FEL  wave  equation,  assume  a  circularly  polarized  plane-wave 
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is  established  in  the  undulator  as  in  2-3  of  the  previous  section  with  a  vector  potential  of  the  form 

E 

A  -  — (sin\|f,cos»|i,0)  ,  2-23 

where  all  the  variables  are  the  same  as  those  defined  in  2-3.  Using  the  vector  potential  in  the  wave 
equation  [Ref.  8] 


r 


2-24 


where  is  the  transverse  ciurent  from  the  transverse  motion  of  the  electron  beam,  and  the  SVAP 
approximation  to  eliminate  all  terms  with  two  derivatives,  gives 


dt^j 


2  dEg 

— :^(cos»|i,-sln\ir,0)  - 
c  dt 


— 5~(sin\|J,cosi|»,0) 
c  dt 


2-25 


The  sin  and  cos  functions  of  y  in  equation  2-25  rotate  fast  with  respect  to  the  slowly 
evolving  complex  envelope.  To  eliminate  this  fast  rotation  in  equation  2-25,  two  orthogonal  unit 
vectors  can  be  defined, 

£,  H  (cosy, -sin\j/,0)  ,  Cj  ■  (sin y, cosy ,0)  ,  2-26 

to  project  into  2-25.  This  reduces  Maxwell’s  second-order  vector  equation,  2-25,  into  two  first-order 
scaler  equations 


V 


dt^J  c  dt 


2-27 
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and 


e,- 


V  -  ±^] 


^  c  at  ' 


2-28 


The  single-particle  current  is  given  by  J^  =  -ecp6’(jf  -  fj)  where  f.  is  the  position  of  the 
i’th  electron,  and  6^(jr-  /p  is  a  three  dimensional  Dirac  delta  function  [Ref.  9J.  This  implies  that  the 
transverse  single-particle  current  with  the  transverse  motion  frpm  2-7  of  the  previous  section 
is 


=  -^-^(cos(A:^),sin(*  z),0)6’(x  -  fp  ,  2-29 

Y/no)^ 

where  B  is  the  field  strength  of  the  imdulator’s  magnetic  field,  m  is  the  electron  mass  and  (0^  is  the 
radial  frequency  of  the  undulator's  magnetic  field.  Substituting  the  single-particle  transverse  current, 
2-29,  into  the  first-order  Maxwell  equations,  2-27  and  2-28,  and  taking  the  appropriate  dot  products 
gives 


ymw. 


cos(^  +  k^z)^\x  -  rp 


2-30 


and 


/,  -^2  =  ~  ^-sin(ti>  +  •  2-31 

Ywo), 

Summing  over  all  the  particles  to  obtain  the  total  transverse  current  and  averaging  both  equations  at 
a  fixed  time  over  a  volume  element  dv  that  is  much  smaller  than  the  coherence  volume  but  much 
larger  than  the  optical  wavelength,  gives  [Ref.  8J 
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l3£, 

c  dt 


2-32 


and 


=  -^^f£p(cos(C  +<!»))  , 
y 


c  dt 


2neK 


p(sin(C  +  <t>)) 


2-33 


where  p  is  the  particle  density  and  the  brackets  represent  the  ensemble  average  over  all  the 
electrons.  Equations  2-32  and  2-33  can  be  simplified  further  by  introducing  another  dimensionless 
parameter 


j  ^  SNjenKLfp  2-34 

y^mc^ 

representing  the  dimensionless  current  of  the  electron  beam  [Ref.  8].  The  dimensionless  current,  j  , 
along  with  the  dimensionless  field  strength  \a\  simplifies  equation  2-32  and  2-33  to 

|al  =  -y(cos(^  +<[))>  , 
and 

^  =  -j^(sin(C  +<!>)). 


The  above  two  equations  describe  how  the  complex  optical  field  evolves  with  the  electron 
beam.  With  high-current  densities,  j  »n,  the  gain  will  be  high  since  |o|  will  be  large.  Also, 
strong  optical  fields  will  slow  down  the  evolution  of  the  optical  phase.  The  above  wave  equation 
used  with  the  pendulum  equation  of  section  B  describe  the  simple  PEL  optical  and  electron  evolution. 
In  general,  the  pendulum  and  wave  equations  have  no  simple  analytical  solution,  but  are  easy  to 
simulate  numerically  to  show  the  major  features  of  PEL  physics.  Section  E  will  present  numerical 
simulations  that  demonstrate  important  PEL  physics. 
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D.  FREE  ELECTRON  LASER  GAIN 


1.  Low-Current,  Weak-Field  Gain 

In  weak-fields  and  low-current,  an  approximate  solution  to  the  pendulum  equation,  2-21,  can 
be  obtained  by  assuming  that  <)>  and  |a|  are  constant.  Using  conservation  of  energy  in  the 
undulator  and  relating  Av  ,  from  the  pendulum  equation’s  solution,  to  the  change  of  energy  in  the 
optical  field,  an  analytic  equation  can  be  found  to  give  the  gain  of  the  FEL  for  weak-fields  and  low- 
current  [Ref.  8]. 

The  derivation  begins  by  finding  a  second-order  perturbation  solution  to  the  pendulum 
equation  assuming  weak-fields  and  low-current  so  that  ^  »  0  and  |a  |  «  0  .  The  solution  is 


C(^)  = 


l^[cos(^j,  +  VpT)  -  cos(Co)  +  v„''S‘n(Co)J+- ••  • 

Vo 


2-37 


so  that 


v(x)  =  Vg  +  _jl[sin(;o  +  Vpt)  -  sin(Co)|  + 


Jr 


v’ 


-l(cos(2Co  +  2Vot )  -  cos(2C„))  +  cos(vj ) 

4 


2-38 


l-Votsin(;^cos(C„  +  VoT )]+... 


At  injection,  the  electrons  form  an  essential  continuum  that  are  evenly  spread  in  phase  space. 
As  the  electrons  evolve  along  the  orbits  in  phase-space  some  will  gain  energy  and  move  ahead  in 
phase  space  while  others  will  lose  energy  and  fall  behind.  This  bunching  occurs  on  the  scale  of  an 
optical  wavelength.  Since  the  opcicai  wave  evolves  slowly  in  tune  and  space,  each  adjacent  2n 
section  of  phase-space  is  essentially  the  same  as  its  neighbors  over  many  optical  wave  lengths 
allowing  equation  2-38  to  be  spatially  averaged  over  a  single  2n  length  giving 
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2-39 


^v)  =  di;^  =  J^rcos(v„x)  -  1  +  IvoTsinCVot) 

Vo  L 

From  the  resonance  condition,  it  can  be  shown  that  Av  =  AnNAy/y  [Ref.  8].  Therefore, 
the  change  of  energy  in  the  optical  wave  can  be  foimd  from  the  average  change  of  energy  of  the 
electron  beam  due  to  conservation  of  energy.  Namely, 


dP  = 


-(pFdV)Ynic^((v>  -  Vp) 
4nN 


2-40 


where  p  is  the  real  current  density,  F  is  the  ratio  of  the  "effective"  cross-sectional  area  of  the 
electron  beam  over  the  optical  beam  and  is  known  as  the  filling  factor,  and  dV  is  a.  microscopic 
volume  element  one  optical  wavelength  long  and  with  the  thickness  of  the  optical  beam. 

The  optical  power  in  a  similar  microscopic  volume  with  a  length  of  one  optical  wavelength 
and  cross  sectional  area  of  the  optical  beam  is  P  ■=  2E^dV I  .  Therefore,  since  gain  is  defined 
as  G  s  dPjP  ,  the  gain  is 


^  _  -pF</V{ywc*((v>-v„)/4;tN]  2^^ 

I2£^dv/8n] 

Substituting  equation  2-39  for  (v)  and  the  dimensionless  current  density  defined  in  2-34  into  2-41 
gives 


j[2  -  2  cos(Vot)  -  Vj,xsin(VoT)] 


2-42 


Figure  2-2  is  a  plot  of  the  low-gain  curve.  It  is  interesting  to  note  tliat  the  optimal  gain 
occurs  at  «=  2.6  rather  than  on-resonance  or  far  off-resonance.  This  is  because  near  resonance, 
the  electrons  can  not  drop  a  large  distance  through  phase-space.  Far  above  resonance,  >  rt  the 
gain  will  oscillate  many  times  along  the  undulator  resulting  in  small  final  gain.  In  this  low-current, 
weak-field  regime,  gain  is  not  exponential. 
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Figure  2-2:  The  weak  field  gain  spectrum  for  low  current. 

2.  High-Current,  Strong-Field  Cain 

A  gain  spectrum,  Gia^y^)  ,  is  easily  obtained  numerically  by  calculating  the  gain  at  the 
end  of  a  simulation  at  each  point  (v^.a^)  with  current  j  .  Figure  2-3  is  the  result  of  such  a 
calculation. 

For  weak-fields  |  S  n  the  profile  in  Figure  2-3  is  essentially  that  derived  analytically 
in  section  1  and  plotted  in  Figure  2-2;  ie  it  is  anti-symmetric  with  a  peak  gain  of  G  =  0.1 35y  at 
Vp  =  2.6  .  How  ever,  as  the  initial  field  strength  increases,  the  peak  gain  of  Figure  2-3  shifts  to 
higher  Vj,  .  This  is  due  to  the  large  initial  separatrix.  whose  height  is  2v/  |n| ,  that  can  now  trap 
more  particles  farther  off-resonance.  The  trapped  electrons  can  then  lose  more  energy  dropping 
through  phase-space.  Also,  the  peak  gain  diminishes  with  higher  field  strengths.  This  is  satiuation. 
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Figure  2-3.  The  gain  spectrum  for  high  current  with  increasing  initial  field  strength. 

and  is  caused  by  the  trapped  particles  in  the  closed  orbits  of  phase-space.  After  trapping,  the  particles 
"slosh"  back  and  forth  in  phase-space  producing  synchrotron  oscillations  as  will  be  described  in  detail 
later  in  section  E.4.  Finally,  note  that  the  gain  spectrum  distorts  and  the  gain  bandwidth  broadens 
in  Vp  [Ref.  8]. 

3.  (rain  Degradation  Due  to  Beam  Quality 

Free  Electron  Laser  design  often  involves  a  trade  off  between  high-current  and  beam  quality 
from  the  accelerator.  Since  the  optical  wavelength  is  usually  set  by  the  application,  the  dimensionless 
current  j  is  determined  by  either  the  real  current  or  the  number  of  undulator  periods.  Increasing  the 
real  current  degrades  beam  quality  by  producing  a  beam  with  more  energy  spread  which  degrades  the 
maximum  theoretical  gain.  Conversely,  increasing  the  number  of  undulator  periods  to  increase 


17 


coupling,  may  decrease  the  maximum  theoretical  gain  since  the  undulator  is  now  more  sensitive  to 
angular  spread  of  the  beam. 


Emittance  is  the  measure  of  the  electron  beam  quality,  and  is  a  constant  given  by  the  product 
of  the  average  radial  and  angular  spread  of  the  electrons  =  y  6,  •  It  can  be  shown  (Ref.  8]  that 
the  average  decrease  of  phase  velocity  due  to  the  emittance  of  the  beam  is 


Av. 


2nN 

TTF 


{K^klyl  ^  . 


2/\2\ 


2-43 


In  order  to  get  an  understanding  how  the  spread  of  phase-velocities  affects  gain,  it  is 
necessary  to  assume  weak-fields  and  combine  the  wave  and  pendulum  equations.  First,  solve  the 
pendulum  equation  for  the  i’th  electron  with  a  first-order  pemurbation  +  v.x  +  ,  and 

substitute  this  solution  into  the  wave  equation.  Integrating  over  all  initial  phase  velocities.  • 

gives  an  integro-differential  equation  governing  the  evolution  of  the  optical  field 

a(x)  =  fVT't'F(TV‘"’’a(x-tO  . 

2  Jo 

where  F(.xO  =  ^dq  j{q)e is  the  characteristic  function  for  the  distribution  j{q)  ,  and  f[q) 
is  the  initial  distribution  of  electron  phase  velocities  v,.  =  about  .  In  the  integral 

equation,  2-44,  all  reference  to  electron  phases  has  been  removed,  but  there  still  exists  the  average 
distribution  of  initial  electron  phase  velocities  [Ref.  8]. 

Two  common  distributions  exist,  the  Gaussian 


/c(^)  = 


exp(-<7/2aG) 

a. 


2-45 


and  the  exponential 


/«(<?)  = 


exp((7/a,) 


for  q<0 


2-46 
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These  common  distributions  have  an  easy  functional  fonn,  but  real  beams  generally  have  additional 
factors  complicating  their  shape  and  their  exact  functional  form  is  unknown.  Still  the  Gaussian  and 
exponential  forms  provide  realistic  insight  to  the  gain  degradation  problem  which  can  be  explored  by 
substituting  one  of  their  functional  forms  into  the  integral  equation  2-44  and  integrating  it  numerically 


ln(l+G) 

0 

G 


Figure  2-4;  The  weak  field  gain  profile  for  low  current  with  an  increasing  Gaussian  energy  spread. 

to  obtain  the  gain  profile  in  Figure  2-4.  Several  features  aie  immediately  evident  from  Figure  2-4. 
First,  the  peak  gain  decreases  steadily  as  the  beam  quality  decreases  with  bigger  a  .  The  peak  gain 
also  shifts  further  off-resonance  with  poorer  beam  quality  and  finally,  the  gain  bandwidth  becomes 
wider.  These  general  features  are  common  to  many  other  beam  distributions,  but  the  details  of  the 
degradation  can  be  dramatically  different.  With  high-currents,  j  >l  ,  and  poor  beam  quality,  tlie 
gain  is  exceeding  difficult  to  predict  and  can  be  off  by  several  orders  of  magnitude  [Ref.  8]. 
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4.  High-Current,  Weak-Field  Gain 


In  the  high-current,  weak-field  regime,  the  integral  equation,  2-44,  can  be  used  to  develop 
analytic  solutions  for  the  optical  field  strength  and  phase.  In  order  to  solve  the  integral  equation, 
consider  an  FEL  with  perfect  beam  quality,  f(x^  =  1  ,  and  starting  on-resonance,  =  0  . 
Taking  successive  derivatives  of  2-44  then  gives,  a  =  ija{x)l2  ,  which  has  a  general  solution 
inv''lving  three  roots  of  the  form  a^exp(a^x)  .  The  coefficients  are  determined  by  the  initial 
conditions  a(0)  -  and  a(0)  =  a(0)  =  0  .  The  exponents  are  complex  roots  of  the  cubic 
equation  a’  -  ijfl  =  0  giving  a  solution  which  is  the  sum  of  all  three  roots 


a(x)  =  _^[exp(y/2)''’(/  +v/^)x/2  +  expG/2)''^(i-V^)T/2  +  exp-i(y/2)''^x 


2-47 


During  the  bunching  process,  X  <  x^  »  <.  1  ,  2-47  reduces  to 

a(x)  =  ao(l  +  »7xVl2  +...)  , 


2-48 


implying  that  the  optical  field  strength  remains  esser...auv  constant  while  the  phase  increases 
proportional  to  x^  .  With  high-currents  oidy  the  fastest  growing  root  dominates  after  the  bunch 
is  complete  giving 

la(x)l  -  .^expO/2)'''/3  x/2  ,  tt)(x)  =  072)''’x/2  ,  2-49 

implying  that  the  gain  is  exponential  until  the  onset  of  strong-fields.  With  large  currents,  /  >  1  , 
the  optical  phase  shift  is  large  as  well. 


E.  NUMERICAL  SIMULATION  OF  FREE  ELECTRON  LASER  DYNAMICS 

The  pendulum  and  wave  equations,  2-22,  2-35  and  2-36  are  three  nonlinear  differential 
equations  that  have  no  easy,  general  analytic  solution.  However,  they  are  easy  to  implement 
numerically  and  can  be  used  to  illustrate  major  features  of  FEL  physics. 


20 


1.  Electron  Bundling 

Figure  2-5,  illustrates  electron  bunching  with  a  monoenergetic,  low-current  beam  with  a 
weak  optical  field  present.  The  electrons  are  at  resonance.  v„  =  ?„  =  0  .  The  large  plot  in  Figure 
2-5  shows  the  evolution  of  the  electrons  in  phase-space.  The  phase-space  plot  is  a  plot  of  the 
electron’s  phase,  C.  versus  the  electron’s  phase-velocity  and  reveals  many  of  the  FEL’s  dynamics. 
The  mitial  positions  of  the  electrons  are  shovwi  in  light  grey  and  become  darker  as  they  evolve.  The 
separatrix  separates  orbits  that  repeat  and  close  on  themselves,  closed-orbits,  from  open-orbits  that 
do  not  close  on  themselves.  The  separatrix  shifts  position  with  time  due  to  the  <|)  term  in  the 
pendulum  equation.  The  separatrix  is  plotted  in  this  simulation  as  the  black  line  in  phase-space.  It 
can  be  shown  [Ref.  8]  that  the  separatrix  height  is  2  jo  .  The  remaining  two  smaUer  plots  show 


***  fel  Phase  Space  Evolution  *** 
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Figure  2-5:  Electron  bunching  with  sample  electrons  on  resonance. 
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gain  and  the  evolution  of  the  optical  wave’s  phase  <|>  .  At  time  x  =  the  electrons  enter  the 
undulator  evenly  spaced  in  phase  with  respect  to  the  optical  wave,  and  begin  travelling  down  the 
undulator  slightly  slower  than  the  optical  wave.  The  electrons  whose  phase  n/2  <  ^  +  <|)  ^  3n/2 
lose  energy  and  begin  to  fall  back  in  phase  space  while  those  electrons  with  -nfl  <  ^  +  <|)  <tc/2 
gain  energy  and  accelerate  in  phase  space.  This  bunching,  under  the  right  conditions,  can  trap 
particles  in  the  separatiix  and  leads  to  substantial  gain  and  subsequent  saturation  as  will  be 
demonstrated  in  section  E.3.  In  this  example,  the  gain  is  minimal  since  approximately  half  the 
electrons  gained  and  half  lost  energy.  The  little  gain  that  exists  is  due  to  the  optical  phase  shift,  <(i  , 
that  causes  the  separatiix  to  shift  slightly  from  its  initial  position  leaving  the  average  position  of  the 
electrons  slightly  greater  than  re/2  . 

2.  Low'Currents 

To  achieve  non-trivial  gain,  it  is  best  to  start  the  electrons  off-resonance  as  illustrated  Figure 
2-6.  Here  the  electrons  were  started  at  v  =  2.6  and  ^  2  with  a  low-current  and  weak  initial 
fields.  More  electrons  were  used  to  show  the  virtual  continuum  of  electrons  present  in  tlie  undulator 
as  well  as  increase  the  accuracy  of  the  simulation.  Note  that  the  most  dense  bunching  occurs  near  ^  re 
where  the  phase  space  paths  are  directed  downward  implying  the  highest  loss  of  energy  for  the 
electrons,  and  henceforth  the  largest  gain.  Also,  note  that  the  electrons  followed  the  open  orbits 
outside  the  separatrix,  and  that  there  was  only  a  minimal  shift  in  the  optical  phase. 

3.  Strong-Fields,  Low-Currents,  and  Saturation 

Figure  2-7  shows  a  simulation  that  demonstrates  saturation  in  a  low-current,  but  strong-field 
PEL  as  might  occur  during  a  single  pass  in  a  PEL  oscillator.  In  this  simulauon  the  electron  beam 
is  no  longer  monoenergetic,  but  has  a  Gaussian  distribution  of  phase-velocities  with  a  width  of 
0^  =  5  .  This  more  realistic  situation  prevents  all  the  particles  from  being  trapped  inside 
separatrix,  which  ultimately  lowers  the  maximum  theoreUcal  gain  that  would  have  been  obtained  with 
a  monoenergetic  beam  as  discussed  in  section  D.l  and  D.3.  The  dots  in  phase  space  represent  the 
final  positions  of  the  electrons  at  the  end  of  the  simulation  rather  than  their  time  evolution  through 
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***  FEL  Phase  Space  Evolution  *** 
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Figure  2-6;  Starting  electrons  off  resonance  increases  gain  for  a  single  pass  through  the  undulator. 

phase-space.  The  electron  dynamics  of  this  simulation  illustrate  several  important  FEL  features.  Due 
to  the  initial  energy  spread  in  the  beam,  many  of  the  electrons  did  not  become  trapped  inside  the 
separatrix  and  are  more  spread  out  in  phase-space  than  those  that  are  bunched  inside  the  separatrix. 
Inside  the  separatrix,  the  electrons  bunch  quickly  due  to  the  strong-fields  and  the  optical  power  grows 
exponentially  as  the  electron  bunch  moves  past  ^+(1)  “  ^  at  about  x  =  1/2  .  It  is  the  peak  to 
peak  height  of  the  separatrix,  4  |a  |  '^  ,  at  this  point  that  determines  the  amoimt  of  energy  that  the 
electrons  can  transfer  to  the  wave  since  Av(t)  =  4  .  This  change  of  phase  velocity  is  enough 

to  move  the  election  bunch  across  the  gain  bandwidth.  The  growth  begins  to  slow  down  until  near  “  ti/2 

where  the  electrons  lost  all  the  energy  possible,  and  flow  into  the  -n/2  <^  +  (|>  <n/2  region  of 
phase  space  where  the  electrons  will  take  energy  from  the  optical  wave  and  cause  the  drop  in  gain. 
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Figure  2-7:  A  low  current  beam  in  strong  fields  at  saturation. 

4.  High-Current,  Strong-Fields  and  Synchrotron  Oscillations 

Figure  2-8  corresponds  to  a  high-current  beam  with  j  =  10000,  the  optical  wave  grows 
exponentially,  as  discussed  in  section  D.4,  until  the  onset  of  saturation.  The  saturation  occurs  with 
similar  dynamics  as  the  example  of  section  E.3.  However,  after  saturation  and  the  downward  trend 
of  tlie  gain  evolution  caused  by  the  electron  bunch  flowing  into  the  -n/2  <  ^  +  <(>  <  n/2  region 
of  phase  space,  the  bunch  continues  to  be  driven  by  the  optical  wave  around  the  closed  orbit  inside 
the  separatrix  and  back  into  the  n/2  <  ^  +  <j>  <  3n/2  region  of  phase  space.  There  the  electrons 
transfer  their  energy  back  to  tlie  optical  wave  which  causes  the  second  increase  of  the  optical  power 
that  is  seen  in  the  gain  evolution.  This  oscillation  in  the  gain  evolution  is  known  as  a  synchrotron 
oscillation.  Also  note  how  rapidly  the  optical  phase  evolves  with  time. 

In  Chapter  lU,  it  will  be  shown  that  the  frequency  of  the  synchrotron  oscillation  is 
=  |o|'^  •  Chapter  III  will  also  develop  a  simple  theory  using  the  concept  of  the  a  single 
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***  fel  Phase  Space  Evolution  *** 
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Figure  2-8:  High  currents  and  strong  fields  cause  synchrotron  oscillations  due  to  the  electron  bunch 
"sloshing"  around  phase  space  following  the  closed  obits  in  the  scparatrix. 

macroparticle  to  represent  the  electron  bunch.  This  will  allow  the  development  of  an  approximate 
solution  of  the  wave  and  pendulum  equations  in  tliis  highly  nonlinear  regime. 

5.  The  Tapered  Undulafor 

This  example  illustrates  the  effect  of  tapering  an  undulator  to  achieve  higher  gain.  At 
injection,  the  electron  beam  had  a  Gaussian  energy  spread  with  a  width  of  ■=  5  ,  and  a  current 
of  j  =  2000.  The  tapering  begins  near  saturation  at  x,  «*  0.6  with  a  tapering  rate  of  5  =  70?: . 
Note  that  most  of  the  electrons  were  trapped  inside  the  scparatrix,  but  that  some  remain  outside  the 
separatrix  and  are  randomized  in  phase-space.  These  randomized  electrons  play  a  minimal  role  in 
the  gain  mechanism  whereas  the  trapped  electrons  bunch  near  the  critical  point  lose  energy,  and  lead 
to  gain.  The  "sloshing"  of  the  electrons  due  to  synchrotron  oscillations  is  still  visible  in  the  gain 
evolution,  but  now  the  closed-orbit  area  of  phase-space  is  relatively  small  due  to  the  tapering  of  the 
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Figure  2-9:  A  strongly  tapered,  high-current  undulator. 

pendulum  equation.  This  causes  the  electrons  to  bunch  near  “  3n/2  so  that  there  is 

continual  growth.  The  electrons  still  "slosh",  but  over  a  much  smaller  area  of  phase-space  which  is 
why  synchrotron  oscillations  play  a  relatively  small  role  in  strongly  tapered  undulators  compared  to 
untapered  undulators. 
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Ill  MACROPARTICLE  MODEL 


A.  INTRODUCTION 

The  high-current,  strong-field  regime  of  the  PEL  is  highly  nonlinear  and  difficult  to  analysis 
quantitatively.  However,  simple  analytic  models  can  be  developed  by  treating  the  particles  trapped 
in  the  closed  orbit  region  of  phase-space  as  a  single  macroparticle.  The  macroparticle  approach  was 
originally  developed  and  used  to  derive  simpler  expressions  of  the  pendulum  and  wave  equations  that 
could  be  used  to  facilitate  PEL  design  [Ref.  10,11].  The  macroparticle  approach  has  also  been 
successfully  used  to  study  the  stability  of  Klystron-type  PEL’s  [Ref.  12],  and  to  explore  the  limits  of 
the  pendulum  model  for  PEL  dynamics  (Ref.  13].  In  this  chapter,  simple  solutions  to  the  wave  and 
pendulum  equations  are  derived  using  the  macroparticle  approach.  The  solutions  allow  quantitative 
analysis  of  the  optical  wave  evolution,  and  provide  useful  insight  into  the  trapped  particle  instability. 

The  macroparticle  model  presented  here  is  first  solved  for  the  tapered  undulator  where 
synchrotron  oscillations  are  small.  Then  a  slightly  more  complex  model  is  developed  leading  to 
solutions  that  describe  synchrotron  oscillations  in  the  untapered  undulator.  Comparisons  of  the  simple 
theory  with  particle  simulations  will  be  interspersed  to  demonstrate  the  validity  of  the  model. 

B.  MACROPARTICLE  MODEL 

In  the  high-gain,  strong-field  regime,  tlie  trapped  electrons  bunch  in  the  closed-orbit  region  of 
phase-space  allowing  them  to  be  treated  as  a  single  macroparticle.  The  untrapped  electrons  will 
randomly  disperse  in  phase-space  and  play  a  greatly  diminished  role  in  the  gain  mechanism.  This 
removes  the  required  averages  from  the  wave  equations  and  simplifies  the  system,  equations  2-22, 
2-35  and  2-36,  from  2n  +  1  equations  to  a  system  of  3  equations, 

X  -  60(t  -Xj)  +  la|cos(^+4>)  , 
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|a|  =  -;Vcos(C +«!►)  . 
|a||=  yr*^(C +«l>)  . 


where  9(x  -  x  j)  is  the  step  function  for  the  tapering  constant,  8  .  The  trapped  particles  are  described 
by  a  reduced  particle  current  jj.  »  0.6y  .  The  macroparticle  remains  trapped  near  resonance,  at  the 
bottom  of  the  potential  well,  so  that  0  . 

1.  Solutions  Without  Synchrotron  Oscillations 

The  near-resonance  assumption  allows  for  the  immediate  solution  of  the  pendulum 
equation,  3- la,  to  give  cos(C+<t>)  =  -8/|c|  .  When  combined  with  the  wave  equations,  3- lb  and 
3-lc  the  result  is 

Jr  5  3-2a 


\a\ 


=J, 


M/l 


1  - 


8^ 


3-2b 


A  properly  designed  undulator  requires  6*  <  laF  to  insure  that  there  will  be  closed-orbits  that  trap 
the  electrons  which  is  consistent  with  the  approximation  in  3-2b.  The  solution  to  3-2a  is  found  to 
be 


3-3 


The  solution  to  3-2b  is  found  by  combining  3-3  with  3-2b  and  integrating  with  the  appropriate  initial 
conditions  yielding 


t|.(x)  =  |aj((l +2>,5(x-x,)/|aJ^)'«-l]/8  . 


3-4 


Without  tapering  the  solutions,  3-4  and  3-5,  simplify  to 


K'c)|'  =  la,  I'  .  + 


3-5 


28 


Figure  3-1  illustrates  a  simulation  of  a  tapered  undulator  with  an  initial  optical  field 
=  10  ,  a  dimensionless  current  j  =  20000  ,  and  a  tapering  rate  of  5  =  270n  .  The 
electrons  are  initiaUy  placed  with  a  random  Gaussian  spread  in  phase  velocity  related  to  an  energy 
spread  in  the  electron  beam  by  4nNAyly  =  Av .  The  random  Gaussian  distribution  about  is  given 
by  /g  (^)  =  exp[-9^/20G  ]/y/2nOg  where  the  i’th  electron’s  phase  velocity  is  v.  =  +  q.  In  Figure 

3-1  the  width  of  the  distribution  is  =  3  .  The  large  left  window  in  the  figure  is  a  plot  the 
electron  phase-space;  electron  phase,  ^ ,  versus  electron  phase  velocity,  v ,  where  each  dot  is  the 
electron’s  final  position.  Note  that  the  trapped  electrons  are  bunched  inside  the  separatrix  while  the 
untrapped  electrons  were  dispersed  randomly.  The  separatrix  separates  the  closed-orbit  region  of 
phase-space  from  the  open-orbit  region  and  is  plotted  as  the  thin  black  line  in  phase-space.  The 
electrons  inside  the  separatrix  are  bunched  near  n  which  causes  the  continual  growth  of  the  power. 
Synchrotron  oscillations  are  visible  as  small  modulations  in  the  gain  evolution  and  are  caused  by  the 
"sloshing"  of  the  electrons  inside  the  separatrix.  The  top-right  window  plots  the  power  of  the  FEL 


Figure  3-1:  High  Current,  strongly  tapered  undulator. 
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over  time  with  the  derived  analytic  solutions  superimposed  in  grey.  The  simulation's  gain  evolution 
is  plotted  in  black.  The  flat  initial  portion  of  the  analytical  solution,  T  »  0  ,  and  the  exponential 
gain  portion,  <  t  <  ,  are  from  the  results  stated  in  Chapter  11,  section  D.4.  The  final  portion 

of  the  analytical  solution  is  the  macroparticle  solution  that  begins  after  saturation  when  tapering 
commences  at  t  »  0.3 .  Sixty  percent  of  the  electrons  were  assumed  trapped  in  this  example.  After 
saturation  when  tapering  commences,  the  power  grows  linearly.  The  bottom-right  window  plots  the 
evolution  of  the  optical  phase  with  the  same  three  regions  of  the  analytical  solution  superimposed  in 
grey.  Note  that  the  phase  grows  linearly  in  time  during  the  exponential  growth  as  predicated.  There 
is  close  agreement  between  the  analytic  solution  and  the  simulation.  A  gain  of  G  =  2  x  Kf’  is 
achieved. 

Figure  3-2  illustrates  a  tapered  undulator  with  a  lower  current  of  j  =  10000  and  a 
smaller  taper  rate  of  5  =  200Jr  ,  but  a  strong  initial  optical  field  =  100  .  The  initial  electron 
beam  has  a  Gaussian  spread  ii  pnase  velocity  with  a  width  of  =  5  .  Saturation  occurs  at 
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Figure  3-2:  High  initial  field  in  a  strongly  tapered  undulator. 
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t  j  <=  0.2  since  the  initial  field  is  near  saturation.  The  gain  grows  linearly  after  saturation.  Figure 
3-3  is  similar  to  the  simulation  of  Figure  3-2,  but  the  initial  field  is  much  weaker,  =  10  .  This 
delays  saturation  until  =  0.4  .  Both  simulations  assume  jj.  =  0.6^.  In  both  figures,  the  simple 
macroparticle  model  compares  extremely  well  with  the  simulations.  Synchrotron  oscillations  are 
visible  in  both  figures  as  the  small  modulation  in  the  gain  and  phase  evolution.  The  analytic  solution 


-n/2  C,  37t/2  0X1 

Figure  3-3;  Weak  initial  field  in  strongly  tapered  undulator. 

can  not  follow  the  synchrotron  oscillations  in  the  gain  evolution  since  the  analytic  solution  assumed 
that  the  macroparticle  was  rigidly  trapped  at  the  bottom  of  the  optical  potential  well.  However,  the 
macroparticle  can  oscillate  slightly  as  is  evident  by  the  closed  orbits  inside  the  separatrix  and  the 
slight  modulation  of  the  gain  evolution.  With  tapered  undulators,  these  oscillations  are  relatively 
small  and  can  be  ignored.  In  untapered  undulators,  synchrotron  oscillations  are  substantially  more 
pronounced  and  must  be  accounted  for.  The  next  section  will  account  for  the  oscillations  of  the 
macroparticle  in  the  potential  well  and  accurately  predict  the  synchrotron  oscillations. 
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Figure  3-4  shows  a  simulation  with  parameters  similar  to  the  LLNL  ELF  experiment,  but 


with  better  beam  quality  [Ref.  14],  The  dimensionless  current  is  j  =  7000  with  a  random  phase- 
velocity  spread  of  width  =  3  .  The  initial  optical  field  is  =  12  ,  and  the  taper  rate  is 
6  =  80n  .  Again,  the  analytic  solution  closely  follows  the  simulation.  The  simulation  and  the 
analytic  solution  both  show  efficiency  of  about  30%  which  compares  well  with  the  experiment. 


***  FEL  Phase  Space  Evolution  *** 
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Figure  3-4:  A  simulation  of  an  FEL  with  parameters  close  to  LLNL  ELF  is  compared  with  the 
macroparticle  model. 

2.  Solutions  with  Synchrotron  Oscillations 

The  macroparticle  model  can  be  extended  to  include  synchrotron  oscillations  for  untapered 
undulators,  5=0,  where  the  effect  of  synchrotron  oscillations  is  more  pronounced.  With  the 
assumption  that  the  macroparticle  is  trapped  at  the  bottom  of  the  optical  potential  well.  ^  =  0  ,  the 
pendulum  equation,  3- la,  can  be  solved  for  an  equilibrium  position 
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3-6 


The  macroparticle  is  now  allowed  to  oscillate  near  the  bottom  of  the  optical  potential  well,  around 
so  that 

C  «  C‘  +  Mx)  =  ^  -  <j>  +  x{x)  ,  3-7 


where  Jt(T)  <  n  .  Assuming  that  the  optical  field  does  not  change  rapidly  after  saturation, 
|a|  ==  ^  ==  0  ,  the  solution  for  x(x)  can  be  found  by  combining  3-7  with  the  pendulum  equation 
3- la.  Using  a  small  angle  approximation  for  x(x)  yields 

^  =  X  =  |nJcos(C* +Jc(x) +(j))  ■=  -|ajx(x)  , 

where  |  was  substituted  for  \a  \  since  \a  |  is  nearly  constant  after  satmation.  Equation  3-8  has 
a  solution  of 


x(x)  =  CoCos(Vj(t  -  T^))  .  3-9 

The  frequency  of  the  synchrotron  oscillation  is  Vj  =  •  The  initial  conditions  for  the 

macroparticle  in  phase-space  arc  determined  phenomenologically  by  comparing  with  simulations.  The 
best  solution  is  found  to  be  «  0.6  and  v,,  =  0  as  expressed  in  3-9. 

Combining  3-7  with  3- lb  yields 

|a|=  -7rCOs(C* +-x(x) +<!))  ^ 

after  the  small  angle  approximation  is  applied.  Combining  3-9  and  3-10  produces  a  solution  for  the 
optical  field  strength 

|n(T)l  =  i!^sin(Vj(x -xp)  +  la^  |  .  3-11 

^  9 
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Combining  the  optical  phase  expression  of  3-2b  with  the  new  optical  field  strength  solution  from  3-11 
gives 


(t>(X)  =  (t)j  + 


^  ^cos(v,(t  -t,)) 

^  VI* 


3-12 


after  expanding  to  first-order  in  x(t)  and  integrating. 

Figures  3-5  and  3-6  compare  the  improved  macroparticle  model  for  untapered  undulators 
with  simulations  using  current  j  =  1000  ,  and  initial  optical  fields  of  |ag|  =  1  and  jOgl  =  30 
respectively.  These  examples  illustrate  the  macroparticle  model’s  excellent  ability  to  predict  the  onset 
of  saturation.  Both  simulations  use  nearly  monoenergetic  beams  described  by  =  1  with  1000 
sample  electrons  and  300  time  steps.  Note  that  in  Figure  3-5,  the  FEL  just  saturates,  whereas  in 
Figure  3-6,  the  FEL  saturates  and  undergoes  nearly  two  synchrotron  oscillations.  This  is  due  to  the 
high  initial  field  that  causes  saturation  to  happen  at  a  much  earlier  time. 
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Figure  3-5:  A  comparison  of  the  macroparticle  model’s  solutions  with  a  weak  initial  field  simulation 
and  current  j  =  1000. 
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Figure  3-6:  A  comparison  with  higher  initial  fields  and  a  current  of  j  =  1000. 

Figures  3-7  and  3-8  demonstrate  the  macroparticle  model  capabilities  with  much  higher 
currents  of;  =  10000,  and  initial  fields  of  (aj  =  10  and  |aj  =  100  .  Both  simulations  have 
nearly  monoenergetic  beams,  0^  =  1 ,  with  1000  sample  particles,  use  600  time  steps,  and  assume 
sixty  percent  of  the  electrons  are  trapped.  In  both  examples,  the  macroparticle  model  predicts 
samration  very  well  and  follows  the  synchrotron  oscillation.  The  analytic  optical  phase  solutions  also 
follow  the  simulation  closely  and  show  the  synchrotron  oscillations.  With  the  higher  initial  field, 
saturation  occurs  sooner. 

Figure  3-9  shows  a  simulation  of  an  untapered  undulator  with  a  much  higher  current  of 
j  =  20000,  and  nearly  perfect  beam  quaUty,  =  1  •  Note  the  very  tight  bunch  inside  the  separatrix 
even  after  nearly  four  synchrotron  oscillations.  Once  again  the  analytic  solution  follows  the 
sunulation  closely,  but  the  predicted  synchrotron  frequency  is  slightly  higher  than  in  the  actual 
simulation.  T ’ate  that  the  analytic  power  solution,  P(t)  ,  after  sannation  is  sinusoidal  but  in  the 
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***  FEL  Phase  Space  Evolution  *** 
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Figure  3-7:  High  current  with  3  synchrotron  oscillations. 

simulation  the  synchrotron  oscillations  are  not  exactly  sinuscndal.  This  is  due  to  the  assumption  dial  |d  |  =  0 
after  satm-ation  which  is  not  completely  accurate.  In  the  simulations,  the  peaks  of  the  power  curve 
are  sharper  than  the  troughs  since  |a(T)|  is  not  completely  constant  as  was  assumed.  In  Figure 
3-10,  Coulomb  affects  have  been  examined  for  this  example  and  appear  to  have  little  effect  on  the 
theory. 

As  a  final  example.  Figure  3-11  shows  the  comparison  between  the  solutions  of  3-7  and 
3-8  and  a  simulation  using  5000  electrons  that  has  parameters  similar  to  the  LLNL  ELF  untapered 
undulator  [Ref.  15].  The  initial  field  strength  is  jopl  =  10  ,  the  dimensionless  current  is 

j  =  3350  and  the  random  phase-velocity  spread  is  =  3  .  From  the  power  evolution,  P(t)  , 
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***  FEL  Phase  Space  Evolution  *** 
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Figure  3-8:  Higher  initial  field  causes  saturation  to  occur  sooner. 


it  is  apparent  that  the  electrons  undergo  about  two  synchrotron  oscillations,  and  that  the  analytic 
solution  closely  follows  the  simulation.  The  analytic  phase  solution  also  shows  synchrotron 
oscillations  as  it  follows  the  simulation. 

C.  CONCLUSION 

The  macroparticle  model  provides  simple  analytic  solutions  to  FEL  dynamics  in  the  highly  non¬ 
linear  regime  with  high-currents  and  strong-fields.  The  results  could  provide  the  necessary  tools  to 
lead  to  multimode  theoretical  studies  of  the  trapped-particle  instability,  sideband  development  and 
chaos  in  the  FEL.  Such  simple  analytic  results  can  be  used  in  the  design  of  futiue  high-power,  high- 
gain  FEL’s. 
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***  FEL  Phase  Space  Evolution  *** 
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Figure  3-9:  A  high  current  example,  j  =  20000,  with  four  synchrotron  oscillations. 
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Figure  3-10:  A  high  current  FEL  with  many  synchrotron  oscillations  compared  to  the  macroparticlc 
model  with  Coulomb  affects  considered. 
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IV  TWO  MODE  THEORY 

A.  INTRODUCTION 

In  the  PEL  oscillator,  the  optical  evolution  can  develop  sidebands  and  become  highly  chaotic 
over  many  low  gain  passes  just  as  the  PEL  amplifier  can  develop  sidebands  and  chaotic  behavior  over 
one  high  gain  pass.  This  highly  non-linear  behavior  in  the  PEL  oscillator  is  difficult  to  analyze  due 
to  the  complexity  that  arises  from  the  multiple  longitudinal  optical  modes  interacting  with  the 
electrons.  One  of  the  consequences  of  this  nonlinear  behavior  is  the  trapped  particle  instability  which 
can  be  analyzed  through  longitudinal,  multi-mode  analysis  [Ref.  16-21].  In  this  chapter,  a  simpler 
model  will  be  developed  that  limits  the  number  of  modes  in  the  oscillator  to  two,  which  is  sufficiently 
simple  to  study  the  stability  of  the  PEL  oscillator,  but  not  excessively  simple  so  as  to  be  unrealistic. 
The  two-mode  model  presented  here  will  provide  a  valuable  tool  to  study  long  term  stability  of  the 
optical  wave,  the  development  of  sidebands  and  the  onset  of  chaos  in  a  PEL  oscillator. 

B.  ELECTRON  AND  OPTICAL  DYNAMICS 

To  study  the  electron  and  optical  dynamics  in  the  oscillator  with  two  optical  modes  present, 
equations  jI  i.;otion  for  the  electrons  must  be  developed  and  then  used  to  derive  the  driving  current 
in  the  optical  wave  equation.  This  can  be  accomplished  in  a  manner  analogous  to  the  derivation  of 
the  pendulum  and  wave  equations  of  Chapter  II,  producing  two-mode  pendulum  and  wave  equations. 

These  two-mode  wave  and  pendulum  equations  will  then  be  used  in  computer  simulations  to  explore 
stability  in  the  PEL  oscillator. 

1.  Electron  Dynamics  and  the  Two-niode  Pendulum  Equation 

The  derivation  of  the  two-mode  pendulum  equation  is  very  similar  to  the  single  mode 
pendulum  equation  of  Chapter  II,  so  only  a  sketch  of  the  two-mode  pendulum  equation  derivation  will 
be  presented.  With  two  optical  modes  present,  the  electron  equations  of  motion  are 
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4-la 


me 


dy  ^ 

dt  me 


4- lb 


where  E  and  B  are  the  electric  and  magnetic  fields  of  the  primary  optical  field,  and  and  B,  are 
the  electric  and  magnetic  fields  of  the  secondary  optical  field.  The  charge  and  mass  of  the  electron 
are  denoted  by  e  and  m,  and  c  is  the  speed  of  light.  The  undulator’s  magnetic  field  is  denoted  by 
S  .  Assume  that  the  primary  optical  field  is  of  the  form 


E  =  £(cosi|t,-sin4r,0)  ,  S  =  £(sini)r,cosi|>,0)  , 


4-2 


and  the  secondary  field  is  of  the  form 


E,  =  £,(cosi|f^,-sun|»,,0)  ,  B,  =  B(sint|i,,cosi|f^,0)  , 


4-3 


where  i|»  ~  kz  -  u>t  +  ^  and  i|ij  =  ik^  -  u>jt  +  <l>5.  *  kg  are  the  wave  numbers,  u)  and  are  the 
frequencies,  and  4>  and  are  the  phases  of  their  respective  optical  waves.  Also  assume  the 
undulator  field  is  helical  of  the  form 


=  B„(cos(AoZ),sin(*pZ),0) 


4-4 


where  is  the  wave  number  of  the  undulator’s  magnetic  field  and  z  is  the  longitudinal  coordinate 
for  the  long  axis  of  the  undulator. 

Substitute  4-2,  4-3  and  4-4  into  the  equations  of  motion,  4- 1 ,  and  solve  for  the  perpendicular 
motion  as  in  the  derivation  of  the  pendulum  equation  in  Chapter  II.  Use  of  the  approximations, 
|£i(l  -  P,)  <  |BJ  and  1£,1(1  -  pp  <  |BJ  (Ref.  8]  gives 

=  -^(cos(*5^),sin(*^),0)  ,  ^5 

where  K  s  <B^/2nmc*  is  the  undulator  ptnameter.  'The  energy  equation  is 
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4-6 


^  =  ——{Ecosik^  +  ♦)  +  E^cosik^  + 
<ft  me  Y 


Now  a  relation  between  the  electron’s  phase  with  re^ct  to  the  optical  wave  and  y  needs 
to  be  derived.  As  in  the  single  mode  pendulum  equation,  start  with  ^  1  ~  Pj.  ~  Pz  ■  Using  the 
perpendicular  motion  solved  for  in  4-5,  solve  for  yielding 


p  »  1  - 


4-7 


2y" 


since  the  electrons  are  relativistic.  Taking  the  time  derivative  yields 

rl 


p  = 


4-8 


Now  define  the  electron  phase  with  respect  to  the  primary  optical  wave  so  that  C  a  (Jt  k^z  -  ur. 
Taking  the  second  time  derivative  of  (  yields 


I 


(*  +  *o)c 


4-9 


when  solved  for  p 


Equating  4-8  and  4-9,  and  using  the  resonance  condition  =  y*/(l  +  K^)  gives 

yl  k  yl 


Y  = 


(k  +  k^c  2k^  2k^c 


4-10 


since  k^<k.  Equating  4-6  and  4-10,  combining  with  the  definition  of  C  and  introducing  the 
dimensionless  time  t  *  ct/L  and  dimensionless  optical  fields  \a]  a  4nNKeLEIy^mc^  and 
a  AnNeKLEjy^mc^  gives  the  two-mode  pendulum  equation 

V  s  C*  =  |a|cos(C  +4>)  +  |fl,|cos(*^  +  i|f^  ^ 

As  before,  "®  "  denotes  djdx  .  The  z  dependence  in  the  pendulum  equation  can  be  eliminated  by 
using  the  approximation  it,  and  the  definition Av^  *  ((k  -  k^/k)2nN  =  ((A.^  -  X)lk^2nN 

where  N  is  the  number  of  undulator  periods.  This  gives 
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4-12 


C  =  |a|cos(C  +4>)  +  |ajlcos(C  +  Av^t  +4>5)  . 

Normally  Av^  is  between  2n  and  lOn  since  the  areas  of  interest  are  -  k)2nNlk  »  a'/*  which 

is  the  typical  range  of  synchrotron  oscillations.  It  is  interesting  to  note  that  by  defining  C  with 

respect  to  the  carrier  wave  (X,(i)) ,  the  symmetry  of  the  two-mode  equation  is  broken. 

2.  Optical  Wave  Evolution  and  the  Two-Mode  Wave  Equation 

The  two-mode  wave  equation  is  derived  in  the  same  manner  as  the  single,  mode  wave 

equation.  The  optical  fields  of  4-2  and  4-3  have  the  vector  potential 

-  £  E 

A  =  ---(sin\lr,cosi|i,0)  +  — (sinili  ,cosi|r„0)  4-13 

k  k. 


Substituting  4-13  into  the  general  wave  equation  2-24  and  applying  the  slowly  varying  amplitude  and 
phase  approximations  (SVAP)  to  eliminate  all  terms  with  two  derivatives  yields  the  left  hand  side  of 
the  wave  equation 


li  =— E(cos»|»,-sint|r,0) 

2  1 

- £4>(sint|r,CO5\|r,0) 

c 

+  -£j(cost|i5,-sinijt  ,0) 
c 

2 

-  -£5<j)j(sint|j,,cosi|r,,0) 
c 


The  rapidly  rotating  ijr  terms  of  4- 14  can  be  eliminated  by  defining  four  polarization  vectors 


=  (cosv,-sin\|/,0)  Cj  =  (siny.cosv.O) 

tsj  =  (cosv,,-sinv,,0)  =  (sin\|fj,cosyj.O) 


4-15 


to  project  into  the  wave  equation  producing  four  equations 
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-(£  +  £jCos(4r  -  tj)  +  -  ♦j))  = 

C  C  '. 

-(-£<i)  +  £jSin(^j  -  -  £j4>sCOs(i|r  - 

c  c 

4-16 

-(£cos(i|r  -  »|»s)  -  £4»sin(4r  -  il»j)  +  £,)  =  -—J^  €„ 
c  c 

--(£sin(4r  -i|r,)-£4)cos(i|»  -  -  £j<|>j))  =  €52 

c  c 

Now  the  right  side  of  the  general  wave  equation  needs  to  be  developed  to  find  the  transva-se  beam 
current.  Since  a  single  particle  has  a  transverse  current,  -  fp  [Ref.  9],  giving 

Jjj  =  (ec£/Y)(cos/c^,sinA^,0)8^’\jf  -  fp  when  combined  with  the  transverse  electron  motion,  4-5. 
Projecting  each  of  the  four  polarization  vectors  into  the  single  particle  current  yields 

'^1  "  ^cos(^  + 

-^2  =  — sin(*^z +\|r)df^>(x -/p 

^  4-17 

^cosikfj^  + 

/ly  '^52  =  “SinfV  + 


Summing  the  single  particle  currents,  4-17,  over  all  the  particles,  equating  with  the  left  hand  side  of 
the  wave  equation,  4-16  and  multiplying  each  of  the  four  resultant  equations  of  motion  with 
AnNeKL^Iy^mc^  yields 

|a(  +  (aj|cos(Av,x  +«|>, -<|>)  -  (a,||/in(Av  x  -4>)  =  -y(cos(C  +  <|>)) 

|a|J  +  |aj|sin(Av,x  +  -  <}))  +  la,  |<f,cos(Av,x  +  <t>,  -  4>)  =  y<sin(;  +  «j))>  ^ 

o  ^  o 

|a|cos(Av,x  +  <j>,  -  «1»)  +  [a |4)sin(Av,x  +({>,-  (|))  +  |oJ  =  -j{cos{k^z  +  <|>,))  ' 

0  ^  0 
-|a|sin(Av,x  +<|»,  -<|i)  +  |a |<j>cos(Av^x  +4>,  -(|))  +  |a,|«|>j=  jisinik^^z  +V,)) 

when  the  dimensionless  time  x  is  introduced  and  j  is  the  dimensionless  current  defm^  aS  before. 
Multiplying  the  second  and  fourth  equatioi^  of  4-18  by  /  and  then  combining  the  fn-st  and  the 
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second,  and  the  third  and  the  fourth  equations  yields  the  two-mode  wave  equations 


since  =  (1  -  AVj/2nA0C('t)  +Av^x  +<j>j.  Consistency  between  the  two  equations  of  4-19 

requires  Av^<  2tcN.  This  eliminates  the  second  of  the  equations  in  4-19  so  that  there  is  only  one 
wave  equation  to  follow 


O 

a 


+  fljC 


iAVjt 


4-20 


3.  Computer  Implementation  of  Two-mode  Theory 

The  two-mode  wave  equation,  4-20,  is  not  immediately  implementable  for  computer 
simulation  since  the  two-modes  need  to  be  separated  into  two  separate  equations.  To  separate  a  fromcr^ 
in  4-20,  one  simplifying  approximation  must  be  made,  the  gain  over  a  single  pass  will  be  small.  This 
limits  the  two-mode  theory  to  PEL  oscillators  with  low  single  pass  gain  since  amplifiers  have  large 
gain  over  one  pass.  This  approximation  allows  the  wave  equation  to  be  integrated  over  time  to 
eliminate  one  of  the  optical  terms.  Before  separating  the  two  optical  terms  in  4-20,  a  value  forAVj 
must  be  established.  Since  the  sideband  is  the  mode  that  is  modulating  the  carrier  wave,  and  is 
produced  by  the  synchrotron  oscillations,  it  will  have  a  frequency  that  is  different  from  tlie  carrier 
wave  by  2jt ,  the  frequency  of  one  synchrotron  oscillation.  Therefore  Av^  «  27t . 

To  separate  the  two-modes  out  of  4-20,  integrate  4-20  over  time  x  from  0  to  1.  using  the 
two  approximations  above.  This  yields 

4-21 

Aa  =  -y(^e"'^<’Wx) 

To  find  the  change  in  the  sideband,  multiply  4-20  by  e  and  integrate  over  time  again  yielding 

An,  =  -/(j['e-«'"-^>dx)  .  4-22 
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Since  it  is  easier  to  deal  with  real  rather  than  complex  numbers  in  computers,  equations  4-21 
and  4-22  can  be  broken  into  real  and  imaginary  parts  by  application  of  Euler’s  theorem  yielding 

{cosQdx 

tsa^  =  yj^  (sin^rfi 

=  -yj^  (cos(C  +27iT)>rfT 
=  yj^  <sin(^  +  27ix)>rfT 


4-23 


Similarly,  the  two-mode  pendulum  equation,  4-12,  can  be  broken  into  its  real  and  imaginary  parts 
^  =  5  +  a^cos^  -  a,sinC  +  cos(^  +2;tT)  -  sin(^  +2jtt)  ,  ^"24 


where  8  is  the  tapering  constant  or  artificial  acceleration  applied  to  the  electrons  after  saturation  in 
order  to  maintain  resonance.  Unlike  an  PEL  amplifier,  the  tapering  constant  is  always  on  in  a  PEL 
oscillator  if  there  is  to  be  tapering  applied. 

The  program  to  implement  the  two-mode  model  is  now  straight  forward.  There  are  three 
main  loops.  The  first,  inner  most,  loop  integrates  the  pendulum  equation,  4-24,  over  all  the  electrons 
and  keeps  a  sum  of  cos^,  and  cos(C,.  +  2nx)  to  be  used  as  the  averages  in  the  wave  equation.  The 
next  outer  loop  is  the  time  loop  repeated  for  each  time  step  of  a  single  pass.  It  integrates  the  wave 
equations  using  the  averages  produced  by  the  inner  most  electron  loop.  The  averages  are  reset  for 
each  time  step.  The  outermost  loop  is  the  pass  loop  that  controls  the  number  of  passes  through  the 
undulator.  At  the  end  of  each  pass,  the  optical  wave  is  attenuated  by  the  oscillator  Q  factor. 

4.  Simple  Oscillator  Model 

Por  low-currents  and  wcak-fields,  the  gain  can  be  calculated  horn  2-42  as  approximately  C  =  0.135y 
for  Vp  =  2.6  for  each  pass  through  the  undulator.  Por  many  passes,  the  power  will  beP  = 


46 


where  n  is  the  number  of  passes.  The  loss  of  optical  power  due  to  cavity  losses  is  related  to  the  "Q" 
factor  of  the  oscillator  so  that  over  many  passes  the  power  in  the  oscillator  will  be 

F  =  F„  exp(0. 135 jn  ~  n/(2)  =  exp (0.1 35/0  -  l)«/0  .  4-25 

From  this  simple  model,  any  combination  of  JQ  =  Cl  and  nlQ  =  C2,  where  Cl  and  C2  are 
constants,  will  scale  the  power  and  phase  evolutions  to  look  identicai  «ttcr  n  passes.  With  this 
simple  model,  many  different  combinations  of  n ,  /  and  Q  can  be  explored  with  one  simulation. 


C.  RESULTS  OF  COMPUTER  SIMULATION  OF  TWO-MODE  MODEL 
1.  Comparison  with  Existing  Simulations 

In  order  to  establish  the  validity  of  the  two-mode  model  and  the  correcmess  of  its  integration 
routine,  simulations  using  the  two-mode  model  over  one  pass  are  compared  with  existing  simulations 
known  to  give  correct  results  over  one  pass.  Figure  4-1  and  4-2  compare  the  two-mode  model  with 
the  single-mode  model  used  to  produce  the  earlier  FEL  simulations.  In  both  figures,  the  large 
window  is  the  phase-space  evolution  of  the  electrons  with  the  final  phase-space  position  plotted  at 
the  end  of  the  undulator  pass.  Over  one  pass,  the  evolution  of  the  second  mode  will  not  have  a 
significant  affect  on  the  power,  phase  or  electron  phase-space  evolution  if  it  is  started  at  zero  power 
since  it  does  not  grow  significantly.  From  Figure  4-1  and  4-2,  it  is  apparent  that  the  phase-space 
e\'olution  is  nearly  identical.  The  optical  power  and  phase  in  both  simulations  produce  nearly  the 
same  final  powers  and  phase.  From  this  comparison  it  can  be  seen  that  the  two-mode’s  integration 
technique  is  probably  correct. 

The  two-mode  model  can  be  compared  with  a  much  more  sophisticated  model  to 
demonstrate  that  the  two-mode  model  accurately  predicts  the  onset  of  a  single  sideband.  Figure  4-3 
shows  a  simulation  of  the  two-mode  model  for  /  =  5  and  0=7  over  KXX)  passes.  Note  how 
quickly  the  main  mode  grows,  saturates,  and  has  a  decaying  oscillation.  The  sideband  has  a  similar 
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Figure  4-1:  Single-mode  simulation  for  comparison  with  the  two-mode  model. 


pattern  but  oscillates  about  a  much  lower  power,  about  1/5  the  power  of  the  main  mode.  Figure  4-4 
is  a  multi-mode  simulation  where  each  mode  is  freely  aUowed  to  evolve  at  its  natural  frequency.  This 
simulation  has  a  random  fluctuation  with  a  rms  value  of  6a  =  5  added  to  the  initial  optical  Held 
strength  to  simulate  noise  in  the  optical  wave.  The  left-center  window  displays  the  evolution  of  the 
optical  power,  la(z,n)l ,  over  1000  passes  through  a  gray  scale  where  black  is  la(r)l  =  0  and  white 
represents  la(2)l  =  25.5  with  a  contour  line  at  la(z)l  =  12.8.  The  window  has  a  width  of  two 
slippage  distances  with  periodic  boundary  conditions.  The  slippage  distance  is  the  interaction  distance 
that  the  electrons  can  couple  with  the  optical  wave.  The  slippage  distance  is  the  distance  that  the 
optical  wave  gains  over  the  electron  pulse  while  overtaking  it  during  the  trip  down  the  undulator  and 
is  given  by,  NX,  where  N  is  the  number  of  undulator  periods  and  X  is  the  wavelength  of  the  main 
lobe  of  the  optical  wave.  For  each  pass,  the  electrons  start  almost  uniformly  spread  in  phase-space 
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Figure  4-2:  Two-mode  model  for  comparison  with  figure  4-1. 


at  each  site  along  2,  but  have  a  small  random  displacement  in  phase  where  a  rms  spread  of 
6^  =  0.0001  is  added  to  simulate  shot  noise.  The  optical  wave  begins  with  no  modulation  but 
undergoes  substantial  gain  during  the  first  100  passes.  The  bottom  left  window  plots  the  gain  per 
pass.  By  n  =  50,  the  optical  wave  is  already  displaying  some  modulation  and  by  n  «  200  the 
modulation  is  easily  visible  as  the  fluctuations  in  intensity  across  the  middle-left  window.  The  top- 
left  window  shows  the  final  power  at  the  end  of  the  1000th  pass  which  is  clearly  modulated.  The 
modulation  is  caused  by  the  electrons  which  slip  back  past  the  optical  field  sites  and  execute 
synchrotron  oscillations.  The  synchrotron  frequency  is  imposed  on  the  optical  wave  over  many 
passes,  so  that  only  a  small  fraction  of  a  synchrotron  oscillation  is  sufficient  to  pass  information  to 
the  optical  wave.  The  bottom-right  window  plots  the  total  average  power  in  the  oscillator.  Note  the 
power  saturates  at  n  =  100,  but  begins  to  grow  slowly  at  n  ~  250  as  the  modulations  in  the  main 
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Figure  4-3:  Two-mode  simulation  with  j  =  5,  Q  =  7,  for  comparison  with  figure  4-4. 


lobe  become  more  pronounced  from  the  growing  sideband.  The  sideband  and  the  main  lobe  are 
clearly  visible  in  the  middle-center  window  where  the  power  versus  frequency  is  plotted  with  the 
same  gray  scale  at  the  end  of  each  pass.  The  main  lobe  is  clearly  visible  as  the  thick,  white  line 
down  the  center  of  the  window  with  the  sideband  approximately  2n  away  becoming  visible  at  about 
n  ~  250.  The  top-center  window  is  a  plot  of  the  final  power  versus  frequency.  The  ratio  of  power 
between  the  sideband  and  the  main  lobe  is  about  1/5  which  compares  favorably  with  the  two-mode 
simulation  of  Figure  4-3.  The  right-middle  window  is  a  plot  of  electron  phase-velocity  at  the  end  of 
each  pass.  It  is  interesting  to  note  the  phase  velocities  spread  out  as  the  sideband  becomes  more 
pronounced.  The  top-right  window  plots  the  final  electron  phase-velocity  distribution.  The  bottom 
center  window  is  a  plot  of  the  gain  spectrum.  In  die  multi-mode  simulation,  the  total  hnal  power  is  3.2  x  1(F 
which  compares  favorably  with  the  total  power  of  3.2  x  10^  in  the  two-mode  simulation. 
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******  fel  Wrap  Evolution  ****** 
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Figure  4-4:  Multi-mode  simulation  with  j  =  5,  Q  =  7,  showing  the  onset  of  a  single  sideband. 
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Figure  4-5  and  Figure  4-6  show  a  two-mode  and  multi-mode  simulation  similar  to  the 
examples  above,  but  now  (2  =  9.  The  higher  Q  value  allows  tlie  sideband  to  fully  develop  to  a 
power  comparable  to  the  main  lobe.  In  Figure  4-5,  the  two-mode  simulation  shows  the  sideband 
being  fully  developed  at  n  =  390.  Before  becoming  fully  developed,  the  sideband  undergoes  two 
large  fluctuations  before  reaching  the  power  of  the  main  lobe  and  stabilizing.  These  fluctuations  will 
be  examined  in  more  detail  in  the  next  subsection.  The  multi-mode  simulation  shows  the  sideband 
developing  much  sooner,  and  being  completely  developed  by  n  «  170  without  the  large  fluctuations 
of  the  two-mode  model.  In  the  multi-mode  simulation,  single-mode  saturation  occurs  at  /;  •»=  50,  but 
growth  of  the  total  power  continues  as  the  sideband  grows  from  n  «  50  to  n  =  170.  A  second 
sideband  begins  to  appear  at  n  -  180,  but  never  grows  to  a  comparable  power.  The  two-mode 
simulation  has  a  total  power  of  P  =  4.4  x  10^  which  compares  well  with  P  -  5.9  x  10^  for  the  multi- 
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***  FEL  Phase  Space  Evolution  *** 
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Figure  4-5:  Two-mode  simulation  showing  sideband  growth  that  is  comparable  to  the  main  lobe 


mode  simulation.  In  the  multi-mode  simulation,  the  optical  power  is  strongly  modulated  from  the 
very  beginning  of  the  simulation  because  of  the  early  growth  of  the  sideband. 

The  last  four  examples  were  not  to  show  that  the  simple  two-mode  model  could  get  the 
same  answer  as  the  much  more  complicated  multi-mode  simulation,  but  were  to  show  that  the  simple 
two-mode  model  is  basically  correct  and  produces  results  that  can  give  insight  into  the  development 
of  sidebands.  In  the  next  subsection,  the  two-mode  model  produces  results  that  give  insight  into  an 
instability  that  leads  to  the  onset  of  sidebands.  Comparisons  with  the  multi-mode  simulation  show 
that  the  instability  is  real  and  not  an  artificial  feature  of  the  two-mode  model.  It  will  be  the  simplicity 
of  the  two-mode  model  that  will  hopefully  allow  greater  insight  into  the  dynamics  of  the  instability 


at  the  onset  of  sidebands. 
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Figure  4-6:  Multi-mode  simulation  with  early  growth  of  the  sideband  to  powers  comparable  to  the 
main  lobe. 


2.  The  Ousel  of  Sidebands 

Figures  4-7  through  4-12  illustrate  a  very  interesting  phenomena  during  the  development 
of  a  sideband  as  Q  is  increased  gradually  over  different  simulations  while  holding  j  =  5  constant. 

In  Figure  4-7,  where  ^  =  5,  the  sideband  only  grows  to  about  1/5  of  the  power  of  the  main  lobe. 

Both  the  main  mode  and  the  sideband  saturate  very  early  at  n  ~  50,  and  the  sideband  has  a  small 
"ripple"  that  quickly  decays  away.  The  main  lobe  achieved  a  power  of  P  =  2.0  x  10^ .  In  Figure  4-8, 
Q  -  6,  and  the  power  in  the  main  lobe  is  slighdy  higher  at  P  =  2.5  x  UF,  but  the  sideband  is  still 
about  1/5  of  the  power  of  die  main  lobe.  The  "ripple"  in  the  sideband  decays  slower  than  in  Figure 
4-7  and  the  main  lobe  is  beginning  to  pick  up  a  similar  "ripple".  The  optical  phase  changes  rapidly. 

The  "ripple"  is  now  very  pronounce  in  both  the  main  lobe  and  the  sideband.  In  Figure  4-9, 
Q  -  7.47.  the  "ripple"  no  longer  appears  to  decay.  In  Figure  4-10,  Q  -  7.5,  the  "ripple"  is  a 
growing  oscillation.  The  mean  power  of  the  sideband  is  still  about  1/5  of  the  mean  power  of  the 
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Figuie  4-8:  Onset  of  the  sideband,  j  =  5,  Q  =  6. 
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main  lobe.  In  Figure  4- 1 1 ,  (2  =  7.8,  several  interesting  changes  happen  in  the  simulation.  First,  the 
sideband's  average  power  grows  and  the  oscillations  in  the  sideband’s  power  evolution  cross  the  main 
lobes’s  power  evolution  twice.  Second,  the  main  lobe’s  average  power  falls.  Third,  the  sideband’s 
power  evolution,  after  crossmg  the  main  lobe’s  power  evolution  twice,  oscillates  to  zero,  then  builds 
up  to  a  steady  power  at  nearly  the  same  power  as  the  main  lobe.  Fourth,  in  all  the  earlier 
simulations,  the  electrons  were  contained  within  a  2n  region  of  phase-space  so  that  ,  and 

only  had  to  cross  the  periodic  boundaries  of  phase-space  once.  In  Figure  4-11,  =>  3n  since  the 

electrons  are  widely  spread  in  phase-space.  Finally,  note  the  difference  in  the  optical  phase  evolution. 
In  the  previous  examples,  the  optical  phase  of  the  sideband  and  the  main  lobe  evolved  rapidly,  but 
at  a  constant  rate  and  together.  In  Figure  4-11,  the  optical  phase  evolves  much  slower  after  the 
sideband  reaches  steady  state.  In  Figure  4-12,  Q  =  8,  the  sideband’s  power  evolution  crosses  the 
main  lobe’s  power  evolution  twice  before  it  reaches  a  steady  state.  Again,  the  optical  phase  evolution 
is  slower  after  the  sideband  reaches  steady  state,  and  the  electrons  have  a  large  spread,  5^  ^  3Jt . 

The  oscillations  and  instability  presented  in  Figures  4-7  through  4-12  are  not  due  to 
numerical  integration  errors.  Numerous  simulations  were  tried  with  many  more  particles  and  time 
steps  and  the  oscillations  still  appeared  with  the  same  shape  and  frequency  and  only  slightly  different 
amplitude.  Near  the  critical  Q  value  that  causes  the  ’’ripple"  to  start  growing  instead  of  decaying, 
the  number  of  particles  and  time  steps  do  affect  the  pass  that  sideband  develops  at.  The  "ripple"  still 
has  the  same  frequency.  The  slight  differences  between  the  simulations  with  more  time  steps  and 
particles  are  most  likely  due  to  numerical  accuracy.  Figures  4-13  and  4-14  show  simulations  with 
200  particles  and  201  time  steps  to  compare  with  Figures  4-3  and  4-11  which  were  run  with  30 
particles  and  30  time  steps. 

The  sequence  of  two-mode  simulations  in  Figures  4-7  through  4-12  can  be  compared  with 
multi-mode  simulations  at  critical  values  of  Q  to  show  that  certain  key  features  exist  between  the  two 
models.  Some  care  must  be  exercised  when  making  direct  comparisons  between  the  two-mode  and 
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multi-mode  model  since  the  two-mode  model  does  not  demonstrate  much  sensitivity  to  noise.  Noise, 
however,  has  a  strong  influence  on  when  the  sidebands  begin  to  develop  in  the  multi-mode  model 
and  sidebands  will  not  develop  without  some  noise  present  in  the  multi-mode  simulations.  Figures 
4-4  and  4-6  of  the  previous  subsection  are  far  enough  away  from  the  critical  values  of  j  and  Q  near 
the  instability,  that  the  noise  sensitivity  does  not  appear  to  affect  the  results  substantially  and  the 
multi-mode  model  accurately  reproduces  comparable  sideband  power  predicted  by  the  two-mode 
model. 

In  Figure  4-15,  the  multi-mode  simulation  is  run  with  j  =  5,  and  Q  =  1.5.  The  simulation 
also  has  a  large  amount  of  shot  noise,  5^  =  0.001  which  tends  to  cause  the  sideband  to  develop 
early.  This  region  of  j  and  Q  is  similar  to  Figure  4-10,  so  that  the  sideband  would  now  have  a 
growing  oscillation.  Since  the  multi-mode  simulation  is  run  to  5000  passes,  there  is  probably  time 
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Figure  4-14:  A  simulation  with  200  particles  and  201  time  steps  at  Q  =  7.8. 

for  the  sideband’s  power  evolution  to  grow  and  overlap  with  the  main  lobe’s  power  oscUlation  as 
happens  dramatically  in  Figure  4-11.  The  large  shot  noise  probably  made  this  interaction  happen  at 
an  earlier  time  tlian  in  the  two-mode  simulation.  This  overlapping  of  the  side  and  main  lobe  power 
evolutions  could  be  the  cause  of  the  sideband  that  appears  to  "grow"  out  of  the  main  lobe  in  the 
center  middle  window  of  Figure  4-15  and  then  "grow"  back  into  the  main  lobe.  This  process  could 
repeat  itself  many  times  as  the  sideband  and  main  lobe’s  oscillations  slowly  grow  and  merge  as  in 
the  two  mode  model.  It  is  also  interesting  to  note  that  the  electron  phase  velocities  became  more 
dispersed  upon  the  first  appearance  of  the  sideband  oscillation  in  Figure  4-15.  This  phenomenon  also 
was  apparent  in  the  two-mode  model.  Note  that  the  total  power  plot  in  the  bottom  right  of  Figure 
4-15  is  constant  after  the  first  1000  passes  implying  that  the  main  lobe  and  the  side  lobe  are 
exchanging  power  as  the  merge  and  re-emerge  with  each  other.  This  will  produce  the  out  of  phase 


59 


******  pel  Wrap  Evolution  ****** 
j=5  6C=0.001  Oq=0  3^=1  V^=2.6 

Q=7.5  6=0  N=100  5a=l 


0  n  5000  -47  v  47  0  n  5000 


Fijuiie  4-15:  The  emergence  of  a  sideband  out  of  the  main  lobe. 


oscillations  that  are  visible  in  the  two-mode  simulation  of  Figure  4-10. 

Figure  4-16  is  a  multi-mode  simulation  with  J  =  5  and  Q  =  7.6  with  no  shot  noise,  but 
a  large  random  fluctuations  in  tlie  initial  optical  power,  5a  =  5.  Here,  the  sideband  starts  growing 
almost  immediately  wliich  is  visible  by  the  slight  grey  trail  that  fonns  early  just  to  the  right  of  the 
main  lobe  in  the  center-middle  picture  of  Figure  4-16.  However,  it  is  not  until  n  ~  280  that  the 
electron  phase  velocity  widens  substantially  which  is  when  the  sideband  is  at  nearly  half  the 
maximum  power  of  the  main  lobe.  The  main  lobe  and  sideband  never  merge  in  tliis  simulation,  but 
they  both  have  oscillations  across  different  frequencies  as  they  evolve.  Note  that  the  total  power  does 
not  become  constant  until  n  ~  400. 

The  final  three  figures  of  this  chapter  demonstrate  tlie  simple  oscillator  model  of  section  B.4 
with  the  two-mode  model.  In  these  examples,  jQ  =  39  and  n/Q  =  128.20513  and  are  constants. 
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Figure  4-16;  Multi-mode  simulation  with  the  sideband  developing  early. 

In  Figure  4-17,  j  =  2,  «  =  2500  and  g  =  19.5.  In  Figure  4-18,  ;  =  1 ,  n  =  5000,  and  2  =  39. 
In  Figure  4-19,  j  =  0.5,  n  =  10000,  and  (7  =78.  Note  all  tluee  simulations  are  similar,  and  the 
differences  are  probably  due  to  numerical  error.  These  figiues  also  are  similar  to  Figure  4-11  where, 
j  =  5,  «  =  1000,  and  Q  =  7.8. 

1).  CONCITISION 

Even  though  no  quantitative  theories  for  the  development  of  sidebands  in  the  FEL  oscillator 
have  been  developed  from  the  two-mode  model,  tire  two-mode  model  was  instrumental  in  the 
discovery  of  the  apparent  instability  that  exists  just  prior  to  the  growth  of  the  fust  sideband.  Together 
with  the  multi-mode  model,  the  two-mode  model  will  provide  a  powerful  tool  for  the  exploration  of 
sideband  development.  It  is  the  simplicity  of  the  two-moue  model  that  will  probably  provide  the 
necessary  intuitive  understanding  of  the  dynamics  of  sideband  development  to  allow  a  simple. 
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Figure  4-17;  Simple  oscillator  comparisons,  j  =  2,  Q  =  19.5  and  n  =  2500. 

quantitative  dieory  to  be  developed  explaining  sideband  development.  The  two-mode  model  can  be 
easily  extended  to  three  or  more  modes  without  excessive  loss  of  simplicity  so  that  studies  can  be 
conducted  oi.  the  onset  of  chaos  as  well  as  the  growth  of  many  sidebands. 
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